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We present an extension of the BELM code (Belmont et al 2008) to investigate the microphysics 
of particle acceleration in black holes accretion disc corona. The updated version of the code 
accounts for the dynamics of resonant slab waves as well as their interaction with both leptons or 
protons. It is found that the proton temperature is an important regulating effect of the stochastic 
particle acceleration process in accretion disk corona. We present a preliminary fit of the high 
soft spectral state of Cygnus X-l. 



"An INTEGRAL view of the high- energy sky (the first 10 years)" 9th INTEGRAL Workshop and celebration 
of the 10th anniversary of the launch, 
October 15-19,2012 

Bibliotheque Nationale de France, Paris, France 



* Speaker. 

' Corresponding author 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



Stochastic acceleration in accreting black holes 



Alexandre Marcowith 



1. Introduction 

Black holes in X-ray binary systems are characterized by two main spectral states (see |jl]]). 
At luminosities close to the Eddington limit the spectrum is dominated by a thermal component at 
keV energies and by a soft non-thermal tail with a photon index T ~ 2.3 — 3.5 extending up to MeV 
energies. In this state the source is bright in soft X-rays and faint in the hard X-rays; it is called 
the High Soft State (hereafter HSS). At low luminosities (about 1% of the Eddington luminosity) 
the soft X-ray component is faint and the spectrum is dominated by a hard power-law tail with a 
photon index T ~ 1.5 — 1.9. In this state the source is faint in soft X-rays and strong in hard X-rays; 
it is called the Low Hard State (hereafter LHS). These spectral states are usually explained by an 
hybrid thermal and non-thermal corona model. The LHS state is usually well fitted by a model in- 
cluding multiple Compton up-scatterings of soft photons by a Maxwellian distribution of electrons 
in a hot (kT ~ 50-100 keV) plasma of Thomson optical depth of the order of unity. The HSS is 
usually interpreted as inverse Compton up-scattering of soft disk photons (UV, soft X) by an hybrid 
thermal/non-thermal distribution in a hot coronal plasma. 

Now the origin of the electron energization is still widely debated. One of the favoured process is 
Coulomb heating by a hot population of ions [||], but magnetic reconnection above the accretion 
disk [||] or turbulent heating by the interplay of magneto-hydrodynamic (MHD) turbulence [||] are 
likely to operate the magnetic field effects are included. This work considers the latter process 
as the main source of electron (and/or positron) acceleration process and investigate the possibil- 
ity to build-up non-thermal particle distribution. We propose a microphysical process describing 
quantitatively the transition between the LHS and HSS spectral states. 

2. Physical modelling 

The hot corona is approximated as a sphere of radius R inside which photons, energetic leptons 
(electrons and possibly positrons) and magnetized plasma fluctuations are in a close interaction. 
We also account for the presence of a thermal population of protons with a fixed temperature 
which will contribute to the electron Coulomb heating. We consider a simplified model in which 
magnetic fluctuations are treated as slab MHD modes propagating forward and backward with 
respect to the mean magnetic field. The waves are present at all wavenumbers as part of a turbulent 
cascade. The turbulent energy is injected at large scales Li„j = R and proceeds towards smaller 
scales through non-linear wave-wave interactions down to scales where the Landau damping by the 
different particle species takes over. Waves are chosen in two circular polarization states left (L) 
and right (R) depending on the clockwise or anti-clockwise electric field rotation direction. Waves 
and particles interact only via gyroresonance (transit-time damping will be accounted elsewhere). 
A wave of frequency CO at a particle wave-number k and a particle interact through gyro-resonance 
described by the following condition [j|] 

co — kv\\=n — , n = ±l, (2.1) 

y 

where £l c and y are the cyclotron frequency and the particle Lorentz factor respectively. Electrons, 
positrons and protons can interact with L and R modes in the lab rest frame. 
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Figure 1: L (red) and R (blue) modes absorption rates (normalized to the electron cyclotron frequency) with 
respect to the normalized wavenumber by a population of thermal electrons (continuous line) and thermal 
protons (dashed line) respectively. Also displayed are the dispersion relations of the two modes in light 
continuous lines. R mode absorption rate by positrons is shown in dotted blue line, the absorption rate of L 
modes is superimposed with the electron curve. 



The wave, photon and particles dynamics are treated by solving three coupled kinetic equations 
(see section ||). The wave absorption coefficient a(k) depends on a kernel describing the wave- 
particle resonance absoiption effect modulated by the particle distribution N(p). The particle ac- 
celeration A(p) rate is similarly calculated with the help of the same kernel modulated by the wave 
distribution W(k). This formalism naturally leads to energy conservation as J dka(k) = JdpA(p). 
To illustrate the respective impacts of waves and particles we first display in figure |T]the absorption 
of L and R waves by a population of thermal electrons of temperature of 100 keV and a population 
of thermal protons of temperature of 1 MeV respectively. As can be seen protons mostly absorb 
L modes whereas electrons (to a lesser degree positrons) and protons compete in the absoiption 
process of R modes. In figure ^| we display the particle acceleration rate by the two wave species 
with respect to the particle kinetic energy . 

3. Numerical modelling 

The waves non-linear dynamics are followed using a diffusion term in a kinetic equation as 
described in [|p. The electron and photon dynamics are followed solving a set of coupled non- 
linear kinetic equations as described in I^J. This work is an updated modeling with respect to [f7|] 
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Figure 2: Electron (continuous line), positron (dotted line) and proton (dashed line) acceleration rate wrt 
the particle's kinetic energy by different wave species (L modes in red, R modes in blue). 



using an accurate calculation of particle acceleration and wave damping processes. We account 
at each timestep for the resonances between the particles and the different wave species (Belmont 
et al 2013, in prep) as described in Eq [2.1| . The code has several inputs: the Thomson optical 
depth T, the magnetic compactness lg = (07- /Snm e c 2 )RB 2 , the disk soft photon compactness l s = 
{Ot /m e c )L S /R, where L s is the disk luminosity, the disk temperature Tbb, the proton temperature 
T p and the turbulence compactness I, = (oj /Snm e c 2 )RdB 2 and the turbulence index q. 

3.1 Case study 

We illustrate the outputs obtained from the updated code using the following fudicial example 
withfl = 510 7 cm, T = 0.5, l s = 51(T 3 , £ B = 51CT 2 , j3 p = 0.1 (plasma beta parameter), = 0.38, 
q = 5/3, l t = 0.025. The figure |3] provides the outcoming particle spectrum obtained after the 
turbulent spectrum has reached the damping scales. One can see a dual electron population (one 
thermal and one supra-thermal). The accelerated electrons do not produced a non-thermal power- 
law because of the important radiative losses that produce a pile-up effect. 

3.2 The proton switch effect 

The protons are very efficient in absorbing L and R modes. It happens that if the proton 
temperature exceeds a certain threshold most of the turbulent energy is used in heating the protons 
rather than accelerating the electrons. This threshold corresponds to a plasma parameter of the 
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Figure 3: Electron distribution function wit particle kinetic energy produced in the fiducial example, 
order of one and can be written for typical numbers adapted for X-ray binaries as: 

T -^^ M ^ XT '' x {wLr s f x {w^)- <31) 

Hence if T p < T t j ir the electron distribution harbours a suprathermal component likewise in figure 
|3[ whereas if T p > T t h r the electron distribution is reduced to a thermal component: no strong 
acceleration does occur. The proton switch effect is likely a natural regulation process in these 
plasmas. 

4. Spectral states of Cygnus X-l 

We provide here a preliminary fit of the source Cygnus X- 1 in the high soft state in figure 
The parameters are the following: x = 0.046, l s = 0.011, l B = 0.19, j8 p = 0.33, l s /l, = 0.2. In that 
case the proton temperature is low enough to permit R modes to cascade down to energies where 
they interact resonantly with electrons. A specific work including fits of Cygnus X- 1 in the LHS 
and the transition between the two states will be presented elsewhere. 

5. Conclusions 

The BELM code [|7|] has been updated to account for the production and the effect of turbu- 
lent magnetic fluctuations over particle acceleration in black hole binary systems. The turbulent 
spectrum consists of slab modes with state of polarization either left-handed or right-handed. Pro- 
tons and leptons can have resonant interactions with each type of modes. It is found that the proton 
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Figure 4: The HSS of Cygnus X-l fitted using the updated BELM code. The reduced x is 231/236. Data 
have been taken from [[| . 

temperature is an important regulating effect of the stochastic particle acceleration process in accre- 
tion disk corona as the entire R mode cascade can get absorbed before being able to interact with 
electrons. In that case the electron distribution is reduced to be thermal. Alternatively if proton 
temperature is small enough hence suprathermal electrons are produced and they can comptonize 
both synchrotron and disk radiation. A preliminary fit shows that the high soft spectral state of 
Cygnus X- 1 can be well reproduced. 

The authors thank A.Zdziarski for providing them with the Cygnus X-l data taken from [Jsj] . 
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